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Summary 

A boundary-value problem governing a three-phase concentric-cylinder model was analytically 
modeled to analyze annular interfacial crack problems with Love’s strain functions in order to find the 
stress intensity factors (SIFs) and strain energy release rates (SERRs) at the tips of an interface crack in a 
nonhomogeneous medium. The complex form of a singular integral equation (SIE) of the second kind 
was formulated using Bessel’s functions in the Fourier domain, and the SIF and total SERR were 
calculated using Jacoby polynomials. For the validity of the SIF equations to be established, the SIE of 
the three-concentric-cylinder model was reduced to the SIE for a two-concentric-cylinder model, and the 
results were compared with the previous results of Erdogan. A preliminary set of parametric studies was 
carried out to show the effect of interphase properties on the SERR. The method presented here provides 
insight about the effect of interphase properties on the crack driving force. 


Introduction 

Advanced composite materials are now being used widely for aerospace applications, and the 
development of reliable analysis methods for these composite structures has contributed to the growing 
use of composite materials. Methods for analyzing the deformation of a composite structure under load 
are well established. However, methods for analyzing strength and the degradation of strength during 
service are less developed. The initial strength of an undamaged composite can be analyzed using 
lamination theory and macroscopic failure criteria applied at the ply level. The effect of constituent 
properties can be included in the analysis using micromechanics. However, it is difficult to account for 
the initial presence of microscopic flaws and the growth of microscopic damage during service. For 
polymer matrix composites, the initial presence and growth of flaws in the region of the fiber-matrix 
interface have a large effect on strength retention. Development of reliable models for the growth of 
damage in this interfacial region would contribute to better understanding of the long-term durability of 
composite structures. Anisotropic and/or nonhomogeneous theories can describe the average macroscopic 
effect of microscopic damage such as fiber breaking, matrix cracking, and interface debonding. However, 
a fracture mechanics approach applied at the microscopic level is needed to account for the initiation and 
growth of microscopic damage. After microscopic damage has grown into a macroscopic defectlike 
interply delamination, failure can be analyzed from a macromechanics level (ref. 1). 

On atomic and molecular scales, damage is treated differently in metals, ceramics, and polymers 
because of the inherent differences in the basic microstructures (ref. 2). In fracture mechanics, the driving 
force for the initiation and growth of damage is the strain energy release rate (SERR). The mathematical 
methods used to calculate the SERR for damage growth in a composite are the same for different types of 
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materials. However, the magnitude of the SERR and the critical SERR for the growth of damage depend 
on material properties and molecular damage processes. 

In the present work, the fracture mechanics problem of a crack in the interfacial region of a composite 
is studied. Although the model could be used for analyzing metal and ceramic matrix composites, in this 
report the model is applied to polymer matrix composites. In a simplified view, a polymer matrix 
composite is composed of two materials (fiber and matrix) with different but uniform properties separated 
by a planar interface. In reality, there is often a region of finite thickness near the interface with properties 
different from the fiber or the matrix. This region is sometimes introduced intentionally. For example, a 
sizing is often applied to fibers to improve weaving characteristics and to improve adhesion between the 
fiber and matrix in the cured composite. Although the sizing material is chemically compatible with the 
matrix material, the mechanical properties of the sizing and matrix are not necessarily equal. Another 
reason for the presence of a third region between the fiber and matrix is the effect of the interface on the 
formation of the matrix material during composite fabrication. For thermoplastic matrix materials, the 
degree of crystallinity near the fiber can be different from that in the bulk region. For thermosetting 
matrix materials, the presence of the fiber surface affects the cure reactions so that the chemical structure 
of the matrix material is different near the fiber. The region of finite thickness between the fiber and 
matrix is called the interphase. The thickness and properties of the interphase depend on the particular 
combination of fiber, sizing, and matrix used as well as the composite processing conditions. Other 
researchers have attempted to model variations in material composition and properties in the interphase 
region of a composite material (refs. 3 and 4). This report does not investigate a specific material system. 
Instead, the fracture mechanics model developed in this report is used for a parametric analysis of the 
effect of interphase thickness and stiffness on the stress intensity factor (SIF) and SERR. 

A fracture mechanics analysis of an interface crack in a three-component system (fiber, matrix, and 
interphase) is performed within the confines of the linear theory of elasticity. The crack problem is of the 
mixed boundary value type and is reduced to a system of singular integral equations (SIEs). Except under 
very special circumstances where a closed-form solution can be found (ref. 5), systems of SIEs are solved 
by a numerical method to obtain the crack-tip SIF and SERR. The SIF and SERR are the forces driving 
damage growth and are related to the effects of applied loads, crack geometry, and material properties. 
The goals of this work are to obtain solutions for crack-tip SIF and SERR and to use these solutions to 
evaluate the effects of material properties and interphase thickness on damage propagation. 

A brief description of the formulation of the problem follows. The first step is to develop some basic 
forms capable of modeling a cylindrical fiber. The cylindrical coordinates r, 0, and v are used. The 
variable v is used instead of the more common z in order to be consistent with reference 5. Results from 
this report reduce to those in reference 5 under special circumstances. In the boundary value problem 
considered, there is a coating between a straight elastic fiber of finite radius and an infinite elastic matrix 
material of different thermomechanical properties. The coating could be the interphase in a polymer 
matrix composite or, more generally, any coating in any composite material. The bonding between the 
fiber and the coating is assumed to be perfect. The bonding between the coating and matrix is perfect, but 
there are axially symmetric cracks on the coating-matrix interface. The applied external loads cause 
traction on the surfaces of the cracks. The goals are to obtain a solution for crack-tip singularities, to find 
the stress distribution in the medium, and to find the energy release rates for the propagation of the crack 
along the interface. 

A special case with only two cracks extending to infinity along the axial direction is shown in 
figure 1 . An infinitely long straight fiber and matrix material extending to infinity both axially and 
radially are joined by a third material called the coating. The bonding between the fiber and the coating is 
perfect. The bonding between the coating and the matrix is perfect in the centered region of contact. In 
figure 1, r is the radial direction, v is the axial direction, and (q l5 v^, (|U t 2 , v 2 ), and (p 3 , v 3 ) are the elastic 
constants of the matrix, fiber, and coating, respectively. There is a geometric symmetry in the medium 
with respect to the x = 0 plane. The fiber radius r f is finite, the radius of the coating r c is finite, and the 
radius of the matrix r m is infinitely large. 

To solve the elasticity problem for the geometry shown in figure 1, first the Green’s function for each 
layer is obtained. Using stress and displacement equations and boundary conditions on the interface of the 
coating and matrix, the integral equations of the problem are derived. These integral equations are 
reduced to a system of SIEs. 
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Uniform far-field stresses applied along the radial direction 


4 r 



I x = 0 

T 


Figure 1 . — Three-concentric-cylinder model with interface annular cracks. L , length of the bonded interface; 
Li , II2, length of the debonded interfaces; Ei , E2, E3, modulus of elasticity for the matrix, fiber, and 
coating; r m , radius of the matrix; ry, radius of the fiber; r c , radius of the coating; v \ , V2, V3, Poisson ratios 
for the matrix, fiber, and coating; pi , p2> P3> elastic constants of the matrix, fiber, and coating. 


Analytical Solution 

Governing Equations 

This investigation deals with the plane strain problem (or the generalized plane stress in which the 
mean state of stress across the thickness is considered, or the quasi-plane state of stress in which the 0 
variable is ignored) for bonded materials. The problem considered here is axisymmetric. There is an 
existing geometric symmetry about the x = 0 plane. The external loads are separated into symmetric and 
antisymmetric components with respect to the x = 0 plane, and the only nonzero stress and displacement 
components in an infinitely long cylindrical domain are given by equations (1) to (6) (refs. 5 to 7). In 
these equations, t is a dummy variable of the Fourier domain, F is the Love strain function, v is the 
Poisson’s ratio, and r is the radial distance. 
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Because of existing boundary conditions and the plane strain assumption with no torsional stresses, it 
is convenient to use the Love strain function F to develop the solution following a procedure similar to 
that in reference 5. The Love strain function F can also be called the Green’s function. In the absence of 
body forces, F is defined by the following biharmonic equation in cylindrical coordinates: 


V 2 V 2 F(r,t) = 0 


where the Laplacian operator in cylindrical coordinates is 


V 2 = 


d 2 1 d 2 
— + — - — t 2 
dr r dr 


The solution of equation (7) is as follows, where p = rt 


( 7 ) 


( 8 ) 


F(r,t) = A^ 0 (p) + J Bp^ 1 (p) + C/ 0 (p)+Dp/ 1 (p) (9) 

The Green’s functions F(r, t) for each layer are 


^i( r , t) = A^o(p) + ^iP^i(p) + c i / o(p) + AP / i(p) 

r c t < p < °° 

(10) 

^(g 0 = ^2^o(p) + ^2P^i(p) _l ~^2A)(p) + ^PA(p) 

0 < p <r f t 

(11) 

^(g 0 = ^3^o(p) + ^3P^l(p) + ^3^o(p) + ^3PA(p) 

r f t < P < r c t 

(12) 


where A u B u C u D Y , A 2 , B 2 , C 2 , D 2 , A 3 , B 3 , C 3 , and D 3 are the unknown constants that are functions of the 
transform variable t. 

So that the solution will be bonded for the first and second layer (matrix and fiber) and because of the 
behavior of the Bessel functions at r = °° and r = 0, some of the constants become zero. For the matrix, 

Ci = 0 and D x = 0. For the fiber, A 2 = 0 and B 2 = 0. Therefore, the solutions for the unknown Green’s 
functions of each layer become 
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F 1 (r,?) = A 1 (?)^ 0 (p) + B 1 (?)p^ 1 (p) 

F 2 (r,t) = C 2 (t)l 0 (p)+D 2 (t)l l (p) 

F 3 (r,O = A 3 ( f )^ o (p) + fi3(0^(p) + C 3(0 / o(p) + ^(0 / i(p) 


(13) 


(14) 

(15) 


The functions F u F 2 , and F 3 , and their derivatives, are substituted into the stress and displacement 
equations (1) to (6), where p = /V, r c < r, < °°, 0 < r 2 < r f , r f <r 3 < r c , r f is the radius of the fiber, and r c is 
the radius of the coating. 

The continuity of the radial stresses (eqs. (1) to (3)), the shear stress equation (eq. (4)), the radial 
displacement equation (eq. (5)), and the axial displacement equation (eq. (6)) at the matrix-coating 
interface and the fiber-coating interface must be applied to find the unknown constants A 3 , B 3 , C 3 , D 3 , A u 
B u A 2 , and B 2 . Using r 3 = r x - r c at the matrix-coating interface and defining a = r c t transforms equations 
(1), (4), (5), and (6) into equations (16) to (23): 


l (r c ,x) = - [ 

n J o 


- A X K 0 (a )■ - (a )■ + (1 - 2Vj )B l KJa ) - (a ) 

a 


cos (xt)dt 


2 r 
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71 


-KJa) KAa) 

a 
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+ Cq 


V 
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J 
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a 


+ Z)q 
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cos (xt)dt 


T rjci (r c , v) = — [ t 3 [-A 1 7f 1 (a)+ 25 1 7f 1 (a)- B l aK 0 (a)- 2v 1 2? 1 7f 1 (a)]sin(x)<i/ L 

71 ^0 


x ra3 (r c ,x) = |j o 

M rl (r e , x) = — — | r [- A, (a ) - B t aK 0 ( a . )] cos (xt)dt 

u r3 (r c ,x) = — — f 2 [-A 3 ^ 1 (a)-B 3 aA r 0 (a)+C 3 / 1 (a)+D 3 oc/ 0 (a)]cos(xO^ 

M ;cl (r c ,x) = A_ f f 2 r i 4 1 ^ 0 (a)+B 1 (a^ 1 (a)-4(l-v 1 )^ 0 (a)lsin(xO^f 
7TLL, •'O 


-A 3 A' 1 (a)+B 3 (2£ 1 (a)-aA' 0 (a)-2v 3 A' 1 (a)) 
+ C 3 I x (a) + D 3 (2(1 — v 3 )/ 1 (oc) + oc/q(oc)) 


sin (xt)dt 




*x3 


I x 1 r 2 [A 3 ^ 0 (a)+B 3 [air 1 (a)-4(l-v 3 )li: 0 (a)] 1 , . 

= ' + c J / 0 (a) +Dj [a/ 1 (a) + 4(l-v J )/ 0 (a)]} s '"(»)‘" 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 

(23) 
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Using r 3 = r 2 = r f at the fiber-coating interface and defining y = r f t transforms equations (1), (4), (5), 
and (6) into equations (24) to (31): 


3 r2 


(r f ,x) = - [ 

7 71 J 0 


- c. 


-^ 0 (y)+- / i(y) 

a 


+ 


D 3 [( 2v 3 - 1 )/ 0 (y ) - y/j (y )]| cos(xt )dt (24) 


J r 3 


(r f ,x) = — f 

7 71 J 0 


n 


-K 0 (y)-^-K l (y) 

1 


+ 5 3 [(1-2v 3 )* 0 (y)-y* 1 (y)] 


+ C, 


-7 0 (Y) + -/i(Y) 


+ D 3 [(2v 3 - 1)/ 0 (Y) - a/, (Y)] 


cos(xt)dt (25) 


T rx2( r /’ X ) = “J 0 { C 2 A(Y) + A [2(! - v 2 )A(Y>+ Y 7 o(Y)]} sin(xO^ 

, 2 rf3 J-A 3 ^ 1 (Y) + B 3 [2/r i (Y)-Y^o(y)-2v 3 ^ 1 (Y)]l . . . . 

rx3(r f’ X) tcJo { +C 3 / 1 (Y)+D 3 [2(l-v 3 )/ 1 (Y)+a/ 0 (Y)]| SmX 

u r2 (r f ,x) = - — [ t 2 \C 2 I l (y)+ D 2 yl 0 (y)]cos(xt)dt 

K(J, 2 Jo 

M r3 (r /; x)=-^-J o r 2 [-A 3 /r 1 (Y)-fi 3 Y^o(y)+ c 3 / i(y)+ Z) 3y / o(y)] cos ( xf )^ 

m x 2 x ) = — J t 2 {<Vo(y)+A[yA(y) + 4 ( ! - v 2 ) 7 o(y)]} sin (xt)dt 
7tfl 2 


( \_2_r 2 J a 3^ 0 (y ) + [y^i ( y ) - 4 (i - v 3 )^ 0 ( y)] 

tch 3 Jo ' | +C 3 / 0 (y) + D 3 [y/ 1 (y) + 4(1-v 3 )/ 0 (y)] 


>sin(xt)dt 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 


Boundary Conditions 

Thermomechanical stresses in the medium without the crack are assumed to be known since they can 
be easily calculated. Hence, by simple superposition, the original problem (fig. 2) can be reduced to one 
in which external loads are first applied to a system with the coating-matrix completely separated (fig. 3) 
and perturbation loads are then applied to close the crack in the bonded region (fig. 4). 

The stress-related boundary conditions at the matrix-coating (region 1 to region 3) interface, where 
r x = r 3 = r c are 


<? r i {r c , x) = O r3 (r c , x) = ph{x - £) ± pd(x + £) (32) 

and 

\xl{ r c >x) = X rx3 { r c >x) = ^{ X ~ y ± 4 § (* + y (33) 
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Uniform far-field stresses applied along the radial direction 
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Uniform far-field stresses® 


;* applied along the radial direction 

Ttm mttmttm mm- 


Matrix (region 1) 


Calculated displacement gap g\ alongSr axis and g 2 along x axis due to far-field stresses 


7 

— 

! Coating (region 3) 

") 

1 

< 

r Fiber (region 2) 

! 

j 

r — • 

1 iS 


▼ 

Figure 3. — Elasticity problem — far-field stress applied with no bonding at the 
coating-matrix interface and a gap equal to the calculated displacements 
g\ along the radial axis and g 2 along the axial axis in figure 2. 


t. 



Figure 4. — Mixed boundary value problem — the unbonded problem of figure 3 
with a displacement perturbation to close the crack in the bonded region. 
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where p and q are the unknown stresses outside the crack. The p and q stresses are zero along the crack 
surface. The displacement-related boundaries at the matrix-coating (regionl to region 3) interface are 

u rl (r c , x ) - u r3 (r c , x ) = ^Sfx- £) ± gfi(x + £) (34) 

and 

u xl (r c ,x)-u x3 (r c ,x) = g 2 8(x - ^) ± g 2 8(x + ^) (35) 


where g x and g 2 are known radial and axial distances between the inner diameter of the matrix cylinder 
and the outer diameter of the coating cylinder while these cylinders are completely separated. 

Both stress and displacement boundary conditions are applied at the fiber-coating (region 2 to 
region 3) interface because there is no debonding. The boundary conditions at the fiber-coating interface, 


where r 2 = r 3 = r f , are given in equations (36) to (39): 



O r3 (r f ,x) 

(36) 

x A r f’ x ) 

(37) 

u,. 3 (r f ,x) 

(38) 

u x3 {r f ,x) 

(39) 


Green’s Function 


There are eight unknown constants and eight equations at the matrix-coating (region 1 to region 3) 
and fiber-coating (region 2 to region 3) interfaces. Substituting equations (16) and (18) into equations (32) 
and (33) establishes the following two equations with two unknowns: 


2 

K 


a 

+ (1 - 2Vj )B l K 0 ( a) - B x aK { (a) 


cos (xt)dt = p(x)8(x-^) 


(40) 


and 


2 p 3 -A l K l (a)+2B l K l (a) 

n Jo -B 1 a^ 0 (a)-2v 1 B 1 ^ 1 ( a ) 


sin(x)df = g(x)8(x - 


(41) 


Making an inverse Fourier cosine transformation to equation (40) results in equation (42): 


-A 1 ^ 0 (a)-^-AT 1 (a) 

a 

+ (1 - 2v, ^^(a) - BflX^a) 


=i;k cos(xt) dx 


(42) 
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Using the properties of the Dirac delta function changes equation (42) into equation (43): 


t 3 -A X K 0 (a) — -K x (a)+(\-2v x )B 1 K 0 (a) - B 1 aK 1 (a) = p(b ) cosfo 

a 


Applying an inverse Fourier sine transformation to equation (41) implies equation (44): 

t 3 \-A l K l (a) + B { aK 0 (a)- 2v 1 5 1 £ 1 (a)] = J g(x)S(x-^)sin(xt)ix 

Because of the properties of the Dirac delta function, equation (44) becomes equation (45): 
t 3 \-A 1 K l (a)+ IBJC^a)- B x aK 0 (a)- 2v 1 B 1 X r 1 (a)] = q(x )sin(^) 

Equations (43) and (45) can be written in matrix form as equation (46) with A } and B l as unknowns: 


(i-2v 1 )^ 0 (a)-a*: 1 (a) AW = 1 3 C ° S ^ 

-ATi(ot) 2 J ST 1 (a)-a/f 0 (a)-2v 1J s: 1 (a)JL jBl ( ? )J ^sin(^) 


Solving equation (46) for A x and B x yields equations (47) and (48) with defined in equation (49): 


A / v _ a f/>fe)cos(^)[ag 0 (a)- 2(l - v^fa)] 

1 r 3 A t + #(£)sin(£f)[(l-2v 1 ).fir 0 (a)-a£ 1 (a)] 

-p(y cos (^)[^i(a)] 
r? ’ A i +9©sin(^) ^o(«)— — ^ 


A! = -a 2 K 0 2 (a) + [a 2 + 2(l - v, )] J fi: 1 2 (a) 


Equations (17), (19), and (24) to (29) are substituted into equations (32), (33), and (36) to (39), 
respectively, in order to find the other unknown constants — -A 3 , B 3 , C 3 , and D 3 — at the matrix-coating 
(region 1 to region 3) interface and C 2 , and D 2 at the fiber-coating (region 2 to region 3) interface. Using 
equation (17) in equation (32) results in the matrix form in equation (50): 


A 3 -K 0 ( a) K x ( a) 


+ B 3 [(l-2v 3 )K 0 (a)-aK 1 (a)]\ 


+ C3 ”A)(°0 + — A(°0 
a 


+ D 3 [(2v 3 -l)/ 0 (a)-a/ 1 (a)] 


cos (xt)dt= 


NASA/TM — 2004-212328 


9 



Applying an inverse Fourier cosine transformation to equation (50) yields equation (51): 


t 3 < 


-KJa) — KA a) 
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+ B 3 [(l-2v 3 )Z: 0 (a)-aZ: 1 (a)] 
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-7 0 (a)+ -7,(00 
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= J p(x)8(x-^)cos(xt)r/x (51) 


Using the properties of the Dirac delta function changes equation (51) into equation (52): 


t 3 { 


-K 0 ( a) — K,( a) 
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+ B 3 [(l-2v 3 )K 0 (a)-aK l (a)] 


+ C. 


-I 0 (a)+-I l (a) 

a 


+ D 3 [( 2v 3 - 1 )/ 0 (a ) ■ - aZ, (a)] 


= p(^)cos(^) 


Using equation (19) in equation (33) gives equation (53): 


-r 

K J 0 
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-A 3 K t (a) 

+ B^K^-aK^a)- 2v 3 ZT 1 (a)] 
+ C 3 / 1 (a) 

+ D 3 [2(l- v 3 )/ 1 (a)+a/ 0 (a)] 


sin(xt )dt = q(x )8(x - 


Applying an inverse Fourier sine transformation to equation (53) yields equation (54): 

^(a)/} = Jo * W 8 (* " ^ M xt ) dx 


f- A 3 K 3 (a )■ + B 3 [ 2K X (a )■ - aK 0 (a )■ - 2v 3 K l (a)]| _ 
| + C3/ 1 (oc) + Z) 3 [2(l — V3 )/j(c)c) + oc/q( 


Because of the properties of the Dirac delta function, equation (54) becomes equation (55): 
3 i-A 3 K 1 (a)+ B 3 [2K x (a)- aK 0 (a)~ 2v 3 K x (a)\\ 


+ C 3 /,(a)+ D 3 [2(l - v 3 )7,(a)+ a/ 0 (a)] J - 


(52) 


(53) 


(54) 


(55) 


Using equations (24) and (25) in equation (36) gives equation (56): 


-K^)--K x (a) 

y 


+ ^ 3 [(l _ 2v 3 )Zf 0 (y)- yZf 1 (y)] + C 3 


-7 0 (y)+-7,(y) 

Y 


+D 3 [(2v 3 -l)/ 0 (y)-a7,(y)] = C 2 


1 


-7 0 (y) + — 7,(y) 
a 


+ D 2 [(2v 3 -l)7 0 (y)-yZ,(y)] (56) 


Using equations (26) and (27) in equation (37) yields equation (57): 

- A 3 K x (y ) + B 3 [2K x (y ) - yK 0 (y ) - 2v 3 £ 1 (y )] + 0,7, (y) 

+D 3 [2(l-v 3 )7 1 (y) + a7 0 (y)] = C 2 7 1 (y)+D 2 [2(l-v 2 )7 1 (y) + y7 0 (y)] (57) 


NASA/TM — 2004-212328 


10 



Using equations (28) and (29) in equation (38) yields equation (58): 


-A 3 K l (y)-B 3 yK 0 (y) + C 3 I l (y) + D 3 yI 0 (y) = C 2 I l (y)+D 2 yI 0 (y) 
Using equations (30) and (31) in equation (39) yields equation (59): 


( 58 ) 


4*o(Y) + B 3 [y^(Y)- 4(1 - v 3 )K 0 (y)] + C 3 / 0 (y) 

+D 3 [y/ 1 (y) + 4(l - v 3 )/ 0 (y)] = C 2 / o(y)+ A[yA(y) + 4 (! " v 2 K(y)] (59) 

From equations (56), (57), (58), and (59), the matrix system in equation (60) is established, with Rl, 

R2 , R3 , and R4 detailed in equations (61) to (64): 


-7. k M 


r (l-2v 3 )K 0 (y) 

-Y^i(Y) 


- 7 o(y) 

+ “ 7 i(y) 


^(2v 3 -1)/o(y) A 

-yA(y) 



^i(y) ' 


27 i(y) ' 

-*iM 

-ya(y) 

7 i(Y) 

+ YA>(y) 


v - 2V 3 ^i(y), 


V - 2 Vi(y), 


Y^o(Y) 

_7 i(y) 

-Y^o(y) 


*oM 


Y*i(y) 

- 4 (1-v 3 K(y) 


7 o(y) 


'Vi(y) 

+ 4(1-v 3 )/ 0 (y) 


Aii 1 ) 


'/a" 

B 3 {t) 


7?2 



7?3 

W4 


7?4 


(60) 


where 


Rl = C, 


-7 0 ( Y) + -A(Y) 
a 


+ D 2 [(2v 3 -l)/ 0 (y)-y/ 1 (y)] 

* 2 = C 2 /j (y) + A [2(1 - v 2 )/, (y) + y/ 0 (y )] 

R3 = C 2 I l (y)+D 2 yI 0 (y) 

R4 = C 2 / 0 (y) + D 2 [y/j (y) + 4(l - v 2 )/ 0 (y)] 


(61) 

(62) 

(63) 

(64) 


The solution of equation (60) yields equations (65) to (68) for the unknown constants A 3 , B 3 , C 3 , and 
D 3 in terms of C 2 and D 2 : 



AY 2 (v 2 -v 3 )(/ 0 (y) 2 -/ 1 (y) 2 ) 

-l + v 3 


(65) 
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( 66 ) 


C 3 (f) — C 2 


* 3(0 = ° 

D if (v 2 - v 3 )[/ 0 (yK(y) + A(yKi(y)] 


-1 + v, 


(67) 


D,(,) 


P 2 (-l + v 2 ) 


-1 + Vo 


( 68 ) 


Equations (25) and (27) are substituted into equations (32) and (33) in order to find the unknown 
constants for the coating layer. The matrix system in equation (69) is established to solve for C 2 and D 2 \ 



iV2~ 

C 2 (i)' 


R5 

_N3 

1V4 

D 2 (t)_ 


R6 


Terms in equation (69) are defined in equations (70) to (75): 

m =- 7 «W+^ 


N 2 = - 


— 1 + Vo 


(-l + v 2 )[(-l + 2v 3 )/ 0 (a)-a/ 1 (a)] 

- Y 2 ( v 2 - V 3 ) / 0 (y) 2 - / 1 (y) 2 ][a^o(a)+ *i(a)] 


a 


N4 = - 


— 1 + Vo 


+ 2 {y 2 (v 2 - v 3 )[a/ 0 (a) - A (a)][/ 0 (y)^ 0 (y) + A ( y)*i (y)]} 
A f 3=A(y) 

(-1 + v 2 )[a/ 0 (a) - 2(-l + v 3 )A (a)] 

+ y 2 (v 2 - v 3 )[-/ 0 (y) 2 + a (y) 2 ]^(«) 

- y 2 ( v 2 - v 3 )a (a)[ a(yK (y) + A (yK (y)] 

R5 = cos (?C)p(C) 
i?6 = 


(69) 


(70) 


(71) 


(72) 

(73) 

(74) 

(75) 


The solutions of equation (69) for the unknown constants C 2 (/ ) and D 2 (t) are shown in equations (76) 
and (77): 
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{\ T ( \ { r \ tr\ [- a/ o( a ) + / i(a)]^(C)sin(fC)] r , . ^ 

|/ 1 (a)cos(?CMC)- k — a ; j[-a(-l + V 3 )j 


(77) 


Since A x (t), B r (t), A 3 (t ), B 3 (t ), C 3 (t ), D 3 (t), C 2 (t ), and D 2 (t) are known, the Green’s functions in equations 
(10) to (12) — -Fi(t), F 2 (t ), and F 3 (t ) — are also known. 


Displacement Equations 

In this investigation, the stability of the interface crack at the matrix-coating interface is examined by 
calculating the SERR at the leading edge of the crack. The fiber-coating interface is considered to be 
bonded. Substituting Green’s functions for the matrix F x (t) and coating F 3 (t) layers into displacement 
equations (20) to (23) yields equations (78) to (81). The variables A x and B x are known from equations 
(38) and (39); C 2 (t ) and D 2 (t) are known from equations (76) and (77); and A 3 (t), B 3 (t ), C 3 (t), and D 3 (t) are 
known from equations (65) to (69). 

u r Jr u x) = — f ^-tcos(xt)dt 

riyi } dr v ; 
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= — f — -t cos(xt)dt 

r3V3 ’ K^Jo dr y ’ 


■ ' [ t 2 -^-cos(xt)dt 

7l(i 3 Jo dp 


^-J 0 ?2 [-^ 3 (0 /i: i(p) _jB 3(OP^o(p) + C 3(0 / l(p) + Z) 3(OP / o(p)] COS ( X 0^ 


Ay 2 (v 2 -v 3 )[/ 0 (y) 2 -/ 


2 . / \2 


J-fr 

7i(l 3 Jo 


— 1 + Vo 


Ci _ Pi y1v 2 -vJ[ayK( y ) + , l ( Y )AT l ( y )]j j|(p) cos[xt)di 


g 2 (-l + V 2 ) 

— 1 + Vo 


pj 0 (p) 


,(r 3 ,x) = — — f (l- 2 v 3 )V 2 +-^—-+-— Fr, sin (xt)dt 

?V3 ’ 7tJLL 3 Jo [ V 3) dr 2 r dr J 3 V ’ 

= — f V 2(l-v 3 )^-5-+2(l-v 3 )-^. + (2v 3 -l)F 3 sin (xt)dt 
rcp 3 Jo [ v iJ dp 2 v 'p dp v 3 7 3 J v ’ 


[A3(*K(p)+B 3 (>)[p« l (p)-4(i-v,)sr.(p)] { . 

-I t < >sin (xt)dt 

Jo +C 3 (0 / o(p) + AW[p / i(p) + 4(1-v 3 )/ o ( P )] 


Ay 2 (v 2 -v 3 ) i 0 {yf-i 


2 , / \2 


— 1 + Vo 


^o(p) 


±-\y + { C2 - si „ w „< 


P2( ,lt. Vi) 1[p / '(p) +4 ( 1 - v »KM] 
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Using the same methods, a solution with a similar form can be obtained for antisymmetric loading 
conditions. The solutions for symmetric and antisymmetric loading conditions can be rewritten as 

equations (82) to (85), where J l n (t,r), j\ 2 (t,r), / 21 (t,r), j\ 2 [t,r), Jh{t,r), / 21 (f,r), / 22 (t,r), 

A l5 and A 3 are defined in equations (86) to (95): 


u r Ar,x) = 


1 


7l(ii 


u r3 (r,x) = - 


ujr,x) = 


1 


7t(I 3 
1 


u x3 (r,x) = 


7t(ii 

1 


7t(l 3 


-^)J 0 J u(t’ r ) 

P£)\ 0 J u(p>) 

P^)J o Jh(t,r) 


cos(t^)cos(xt) 

sin(^)sin(xt) 

^cos(^)cos(xt)^ 

sin(t^)sin(x/ L ) 

^cos(^)sin(xt) ^ 
-sin(/^)cos(xt) y 

cos(/^)sin(xt) 

-sin(/^)cos(xt) 


dt 


- + 9(^)J 0 Ju{^ r ) 


dt 

— + 
t 


dt 


<7(^)J o J n{^ r ) 


j+q{^)\ 0 Y 22 (t,r) 


dt 

— + 
t 




sin(t^) cos(xt) 
-cos(t^) sin(xt) y 

sin(t^) cos(jtf) 
-cos(t^) sin(xt) 

sin(t^) sin(xt) 
cos(t^) cos (xt) 

sin(t^) sin(xt) 
cos(t^) cos(xt) 


dt 

t 

dt 

t 

dt 

t 

dt 

t 


(82) 


(83) 


(84) 


(85) 


where 


7 M^O = T{[“ a2 ^(«) +2 ( 1 - v iKi( a )]^(p) + Ki(«)]p^o(p)i 


( 86 ) 


/ 1 1 2 (f,r) = -J-{[(i_2v 1 )a^ 0 (a)-a 2 ^ 1 (a)]^ 1 (p)+[a/r 0 (a) + /r 1 (a)]p^ 0 (p)} (87) 

^1 L 

• / 2i(^'-) = 7-{[ a2 ^(«)- 2a ( 1 - v iKi( a )]^o(p)-«^i(a)[p^i(p)-4(l-v 1 )^o(p)]} (88) 

4(f^) = 7-{[( 1_2Vl ) a ^o(a)-a 2 ^(a)]^o(p)-Ko(«) + ^i(«)][p^i(p)-4(l-v 1 )^ 0 (p)]} (89) 
(-1 + v 2 )p7 0 (p)/i(cc) + / 1 (p)(l - v 2 ) 


J n( r ’ t )= — 


a/ 0 (a)-2(-l + v 3 )/ 1 (a) 

+ y 2 (v 2 - v 3 )/ 0 (y) 2 - ^(y) 2 )^(a) 
+ Y 2 (v 2 - v 3 )/ 1 (a)[-/ 0 (y) 2 + AM 2 ]*^) 


(90) 
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Integral Equations 

The stress-related boundary conditions at the matrix-coating interface are given by equations (96) 
and (97): 

Grl{ r c’ X ) = G r 3 ( r c’ X )= p{ x ) x e (L + Z/) (96) 

\xi{ r c’ x ) = ^ 3 { r c’ x ) = v{ x ) xe(L + L') (97) 

where L is the length of the bonded surface and Z/ is the length of the crack surface. The functions p(x ) 
and q(x) are zero on the crack surface and unknown on the bonded surface. The displacement-related 
boundary conditions at the matrix-coating (region 1 to region 3) interface are given by equations (98) 
and (99): 

U rl( r c’ X )- U r3( r c’ X )= §l{ X ) XeL (98) 
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u xl {r c ,x)-u xi {r c , x ) = g 2 {x) 


X E L 


(99) 


For the symmetric case, 

p{ x ) = p{-x) 
q(x) = -q(-x) 

gi{x) = g l {-x) 

82{ x ) = -82{- x ) ( 10 °) 

For the antisymmetric case, 

p{ x )=-p(- x ) 

q{x)=q{-x) 

8 ih) = -ft(-^) 

8 2 { x ) = 82 {- x ) ( 101 ) 

The problem is solved if p(x ) and q(x) on L are determined. In order to find these functions along the 
bonded surface L, the integral equations are obtained by substituting equations (81) to (84) into equations 
(99) to (100). For the symmetric loading condition, the integral equations (99) and (100) become 
equations (102) and (103): 

f jj - Zllifld + Zilfld [ C os(^)cos(xf)]— }p(^)^ 

1 Jo [ Jo L Hi ^3 J t J 

+ J 0 |j o /|2 ^ ^ + ' /|2 ^ ^ [sin(^)cos(xO]y|?(^)<^ 

J jj J n {t’ r ) _ J 2i{t’ r ) [cos( ^ )sin(x0 ]^l ( ^ )4 

1 Jo [ Jo h F 3 J t J 

+ lo | Jo “ [sin(4) sin (xt)\ y (£)</£ 

For the antisymmetric loading condition, the integral equations (99) and (100) become equations 
(104) and (105): 

j;{j; + [sin(^)sin(xf)]yj^)4 

+ Jo | Jo [- cos (t%) sin(xt )] yj? (£)<^ 
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dr AM [-sin(£)cos(*0]yW>^ 


lim — 

r^a 71 


■ru 


4M _ 4Ml[ cos( ^) cos(x0 ] A qC ^ 


=g 2 ( x ) 


For the combined symmetric and antisymmetric loading condition, the integral equations are given by 
equations (106) and (107): 


ru 


lim 0^ 

r^>a ZK 


Jlijtp t {njtp cos(4)cos(;tf) dfl 
IJ-! ja 3 + sin(f^) sin(xO t \ 


r x 


j\ 2 (t,r ) Ji 2 {t,r ) sin(f^)cos(xf) 


- cos(f^)sin(xO I f 


= &(*) 


q(^)d^ 


ru 


lim ~Ott 

r—>a ZK 


Aijtp Jh{t,r) cos(^)sin(xO 

IJ-! (J-3 _ - sin(?i;) cos(xO_ t 


J 22 (t,r ) J\ 2 [t,r ) sin(^)sin(xf) 1 dt 


■ru 


+ cos(/^)cos(x?) I t 


= g 2 (x) (107) 




The upper and lower kernels correspond to symmetric and antisymmetric cases, respectively. 
Considering the symmetry properties of the functions p(x), q(x), gi(x), and g 2 (x ) — which are given by 
equations (100) and (101) — and using trigonometric identities enables the integral to be rewritten as 
equations (108) and (109): 


ru 


'im t- 

r^>a Zjl 


J\ x (t,r) Ju(t,r)\ 

— -H — - cos(^-x) — >p(Qd!t> 


f” J f“ J \ 2 d\ 2 {t,r^ dt\ „ c 

■ M — — + — - — L sin(^-x) — 

Jo Jo m |x 3 t I 


= &(*) 


r | r r^_^i *| 


lim ~Ott 

r—>a ZK 


4M _ Mid j[ cos( ^ _ x) ] x*l (M 


= £2 00 
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Equations (108) and (109) become equations (110) and (111) with the terms defined in equation 

( 112 ): 


lim 

r—>a 


2n 


J o U" o tf n (/,r)cos(£-x)/y pCQdl, 

+ J Q J o H n (t,r)sm(^-x)t— 


= gi(x) 


( 110 ) 


and 


lim 

r—>a 


2k\ 


“Jo {Jo H 2 i {t,r)[sm(^-x)]—jp(^)d^ 

+ Jo j Jo H 22 {^ r )[ cos (^>- x )]— |tf(£)^ 


= g 2 ( x ) 


(111) 


where 


H n (t,r) 

H n (t,r) 

H 2l (t,r) 

H 22 (t,r) 


J' u (t,r) | Ju(t,r) 

(Xi | i 3 

! ^12 (^’ r ) 

^2l(^’ r ) | ^2l( f ’ r ) 

Hi (I 3 

/ 2 2 (f,r) / 2 3 2 (l,r) 

1^1 H 3 


( 112 ) 


Since the inner integral for t cannot be integrated using usual integration methods, an asymptotic 
approach is applied. When the asymptotic expansions of H^t, r) are considered as t — > 0, some of the 

kernels will diverge. For divergent kernels to be avoided, the Ht factor must be removed. For this reason, 
equations (110) and (1 1 1) are differentiated with respect to x. Using the identities in equations (113) and 
(114) changes equations (111) and (112) into equations (115) and (116): 


— cos(a - x)t = t sin(a - x)t 
dx 


lim 

r—>a 


2n 


— sin(a - x)t = —t cos(a - x)t 
dx 


[ [ H n (t,r)sin(£-x)dt ptfydl, 

JO Jo 

“Jo [jo ^12 (l’ r ) COS (^ <?(^)d£, 


= g'(x) 


(113) 

(114) 


(115) 
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= g' 2 (x) 


(116) 


lim 

r—>ci 


2n 


Jo 1 Jo 
+ 


J <jj // 21 (/ L ,r)[cos(^ - x)]dt'jp(fyd^ 

1 {Jo H 22{t’ r )h in (^- X )] dt \<l(fy d ^ 


Integration of equations (115) and (116) using an asymptotic approach gives equations (117) and 
(118) and the terms in equation (119): 


1 

iim 1 

r->a 271 


Um i 

r->a [ 271 


rW lP £)d^- f W 2 q(0^ 

JO J 0 

f W 3 p(^+ f ~W 4 q(Z>)<% 

Jo J 0 


= &'(*) 


= g' 2 (x) 


(117) 


(118) 


where 


W x = J (^» r )] sin(^ — jc)^ + J* 7/ 1 °j(^r)sin(^- x)dt 

W 2 = \ q [H l2 (t,r)-H~ 2 (t,r)]cos(Z } -x)dt + J H~ 2 (t,r)cosfc-x)dt 
W 3 = J [/7 21 (l,r)-77”(f,r)]cos(^-x)Jf + | o [t, r) cos(L,~ x)dt 

W 4 = (t,rj\sm£-x)dt + (t,r)sin£- x)dt (119) 


Singular Integral Equations and Fredholm Kernels 

In the problem considered here, the bonding stresses p and q are unknown, and g x and g 2 are the 
known displacement differences. Therefore, the displacement-related equations (98) and (99) are used to 
find p and q. From these two boundary conditions, two integral equations are constructed to find the two 
unknowns. In equations (117) and (118), sine and cosine terms are the steady-state oscillating functions. 
Therefore, H^t, r) determines the behavior of the integrals for these equations. For the integrals in 
equations (117) and (1 18) to be calculated for a given interval, the integrand must be integrable over the 
corresponding domain. If the constant terms at infinity in H^t, r) are subtracted from H^t, r), the 
functions H^t, r) - r) become integrable. After the asymptotic values are subtracted and added, the 
integral equations become SIEs with two parts. The dominant parts contain the singularities, and the 
regular parts are the Fredholm kernels. After the asymptotic part is subtracted, the kernels can be 
integrated numerically. But this subtracted quantity has to be added back to the whole equation. That is 
where the dominant part of the SIE comes from. If the limit is taken as 8 — > 0 and r — > r c in the integral 
equations, the Fredholm kernels and the SIEs in equations (117) and (118) become equations (120) and 
( 121 ): 


{-Tf’qM - ^ J* [*i i ( *. WO + j =g,(i) 

-y*q(x)-—f -^-dt - — [ \k n (xj)p(t) + k n (xj)q(t)\lt = Q 1 (x) (120) 

7 \Ja t~X Tl J a 
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773 — t|y*pW--J ^p-dt + -\ [k 2l (x,t)p(t) + k 22 (x,t)q(t)]dt\ = g' 2 (x) 

TC di t X K a 

y * p(x)-—^ q^ dt + — J [k 2l (x,t)p(t)+ k 22 (x,t)q(t)]dt = Q 2 (x) 


where a<x<b and a<t<b. 


y* = 77- {K 1 - 2v i X 1 - v 3 )h - [(1 - 2v 2 )(i - v 3 )]p. 1 } 

M'O 


Ho = 2 ji 3 (1 - v !)( 1- v 3 ) + 2jii (l-v 2 )(l-v 3 ) 


M x ’0 = J o 




{ 2(1 — Vi )oc[ J fiT 1 (ot)] 2 } 


Ho 




A! 


2 ( 1 - Vl ) 


sin[(?-x)s]df 


Ho 


2«(-l + v 2 )(1 - v 3 ) 2 [A(cx)] ^ 


-v 2 )(l-v 3 ) 


kn{ x d) ~ ^2i( x ’0 = 


r 

Jo 


( 1 -V 3 )|4 3 

1 

Ho 

A! 

V 

( 

_ Hi_ 

(1- 

Ho 

A 




ja 2 [£ 0 (a)] 2 + 2(1 - Vj JaATofaJffjfa 
-a 2 [^(a)] 2 

a 2 (- 1 + v 2 )[/i(a)] 2 

-Y 2 (v 2 -v 3 )[/ 0 (y)f 
■< + 2a(-l + v 2 )( 1 -v 3 )/ 0 (a)/i(a) 

-a 2 (-l + v 2 )[/ 0 (a)] 2 
l +Y 2 (v 2 -v 3 )[/ 1 (y)] 2 


■ (l - 2Vj ) 


- (! - 2v 2 )(! - V 3 ) 


yj 


cos[(?-x)s]d,s 


( 121 ) 

( 122 ) 

( 123 ) 

( 124 ) 


( 125 ) 
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Verification of the Singular Integral Equation 

For verification of the analytical result, the SIEs for the three-concentric-cylinder model can be easily 
reduced to the model with two concentric cylinders developed by Ozbek and Erdogan (ref. 5). Making the 
material properties of the coating equal to the material properties of the fiber (i.e., p 3 = q 2 , v 3 = v 2 , and 
E 2 = E 3 ) in equations (120) and (121) yields exactly the same SIEs as in equation 14 of reference 5, as 
shown in equations (130) and (131): 

-y 2 \ dt — — f k n (xj)p(t) + k n (xj)q(t)\lt = Q l (x) (130) 

71 t~X J 

y \ dt + — f \k 21 (xj)p(t) + k 22 (xj)q(t)\it = Q 2 (x) (131) 

where the kernels / 2 are reduced to the one in reference 5 as shown in the following equation: 
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a<x<b a<t<b 


Ho = 2n 2 (1 - v t )(1 - V 2 ) + 2n t (1 - v 2 )(1 - v 2 ) 
Ho=( 1 -v 2 )|io 

yl = ^ Vz ^ [(l - 2v, )(i 2 - (l - 2 v 2 )jiJ 

y *2 = 4- [(l - 2v, )n 2 - (l - 2 v 2 )^i 1 ] 

Ho 


Y 2 =Y 2 


(i-y i yJ{ 2 < 1 ~ v 'H g .( a )] 2 } 


( 1 -V 2 )Ao 


- 2(1 - Vj ) 


( 1 — v 2 )Hi 1 2a(-i + v 2 )(i-v 2 )[/ 1 (g)] 2 


( 1 - v 2 )h 0 I (-i+v 2 )a 2 


- 2(1 - v 2 ) 


k n {x,t) = k n {x,t) 


A,[l + 2 ( 1 -V,^ 0 (a)*r,(«)j UJ J 

L -m^) 

a 2 (-l + v 2 )[/j(a)] 2 

( 1 -V 2 V 1 p7~u; < + 2a(-l + v 2 )(l- v 2 )/ 0 (a)/ 1 (a) 

(!- v 2 )fLo - ot 2 (-! + v 2 )[/ 0 (cx)] 2 


- (l - 2 v 2 ) 


k l2 ( x ,t^) = k 2l [x,t^j = k l2 (x,t) = k 2l (x,t^) 
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(l-v 2 )|X 2 [ 1 (3-2 Vl )a^ 0 2 (a) 

(l - v 2 )A 0 \\[ + 4(1 - v, )K 0 (a)K l (a) - a K\ (a) 


-2(l-v.) 


=f 


+ 


(t-vzK 


(l-v 2 )Ao 


( 1 + V 2 )^2 


a(-l + v 2 )(3 - 2 v 2 )/,) (a) 

-4(-l + v 2 )(l-v 2 )/ 0 (a)/,(a) 

-^i a2 (- 1+v 2K(«)] 


sin[(f- jc),?]^ 


- 2(l- v 2 ) 


k 22 {x,t^) — k 22 [x,t ) 


(136) 


a = r f t y = r f t 

\ = -a 2 K 0 2 (a) + [a 2 + 2(l - v, )]^, 2 (a) 

A, = A, (137) 

A 2 = oc 2 (-l + v 2 )/ 0 (a) 2 -(-l + v 2 )(2 + a 2 -2v 2 )/ 1 (a) 2 

A 2 = (-1 + v 2 )A 2 (138) 

Therefore, the validity of the left side of the SIEs (eqs. (120) and (121)) are established. Reducing the 
right side of the SIEs (eqs. (120) and (121)) to the case of the two-cylinder model (as in eq. (16) of ref. 5) 
yields equations (139) and (140): 


0 = 2(l-v 2 )Mi s;w 

q,=2^ s ;(x) 

a= a (139) 

0 2 =2(1-v 2 )M^( x ) 

Q 2 =2^^g' 2 (x) 

Q 2 =Q 2 ( 140 ) 

The validity of the left and right sides of the SIE is thus established. 
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Numerical Solution 


Fundamental Function of the Singular Integral Equation 

The system of SIEs is derived for three concentric cylinders with two annular interface cracks. In both 
SIEs (eqs. (120) and (121)), the bonding stresses p and q are unknown and g x and g 2 are the known 
displacement differences. For the problem at hand, the SIEs have a simple Cauchy kernel and lengthy 
Fredholm kernels, and it is most convenient to use a numerical method to obtain the solution. Numerical 
solution methods for SIEs have been studied quite extensively (see refs. 8 to 12). There are two methods 
that have been commonly used to solve the integral equations with Cauchy-type singularities. These are 
the quadrature method and the series expansion method. In this study, the series expansion method is 
adopted (see refs. 9, 12, and 13). 

Normalization of the interval of the integrals in the SIEs (eqs. (120) and (121)) yields equations (141) 
and (142): 


-y*q(s) - f J_ ( ~ J_j lA i ' ( 5 > £)/>(£) + k n ( s > = Ci 0 ) (i 4 1) 

Y *p(s) - f f J J [ /c 2 i 0, £)p(£) + k 22 ( s , £)? (£)]cf£ = Q 2 (s ) (142) 

The system of SIEs in equations (141) and (142) can be written as one complex equation. Multiplying 
equation (141) by -i and adding the result to the second equation (142) yields the single SIE in equation 
(143), with the terms defined in equations (144), (145), (146), and (147): 

Y \(s) ~ + K 2 (sX)U^ = QLs) (143) 

<M 5 ) = K y ) + *<?( 5 ) <t>l(^) = p(^) + i#) <t>2 (S) = P(S)-*?(S) ( 144 ) 

Q(s) = Q 2 (s) + iQ l (s) (145) 

p S) = | [k 21 (s£) + k n (s, £)] + i[k 22 (s,q) + k u (5, ^)] (146) 

K P 5 ) = \ [*21 5 ) - k \i £)] - *[*22 Q ~ k \ 1 S)] ( 147 > 


The kernels K ] and K 2 are bounded and (f>, ^), K 2 (s, ^), Q(s) are complex functions. The 

singular behavior of the function (|)(^) at ±1 is determined by the dominant part of the SIE. The function 
(|)(£) satisfies a Holder condition on every closed subinterval of (-1, 1). Its behavior near the ends can be 
represented by equation (148): 


(|)(s) = w(s)/(s) PI < 1 (148) 

The function/is Holder-continuous on any closed subinterval of (-1, 1). This means that p(s ) and q(s ) 
are continuous on the open interval -1 < s < 1 and that (|)(^) has integrable singularities at s ± 1. The 
Fredholm kernels are also Holder-continuous with respect to both variables along the integration interval. 
The dominant part of the SIE is considered to be equal to zero so that the fundamental function can be 
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found. The fundamental function w(s) is obtained from the homogeneous dominant part given by 
equation (149): 

f W ^ 7 d £, + y*w(s) = 0 (149) 

ni J -i s-c, 

The fundamental function w(s) in equation (149) can be expressed as in equation (150) (ref. 13): 

wO?) = (l-sf(l + s) p (150) 

where Xand P are defined in equations (151) and (152): 


X = r co 

2 


P — f i to 

2 


(151) 

(152) 


The term (0 in equations (151) and (152) is defined in equation (153): 


co=ilog| 


r, *\ 

1 + y 

* 

1-7 


(153) 


The boundaries for the real part of Xand Pin equations (151) and (152) are defined in equation (154): 

-1 < Re X < 0 
-1 < P< 0 


(154) 


In equation (150), ir(.v) is the weight function that depends on the index of the singularities. The index 
of the singularity K for a given problem is defined in equation (155): 

K = -(X + p) 


K = — 


\( 1 . ) 


f i . Yl 


+ 


LI 2 J 


l 2 )\ 


K = 1 


(155) 


Solution by Jacobi Polynomials 

After the fundamental function w(s) of the SIE is obtained, the solution of equation (143) becomes 
equation (156): 

^( 5 ) = Xc„w(s)Py p) (s) (156) 

0 


NASA/TM — 2004-212328 


27 



In equation (156), c n are defined as undetermined constants that need to be calculated with 

n- 0, 1, .... The term P^\s) is the Jacobi polynomial. The fundamental function w(s) of the SIE in 
equation (156) is the weight function of the Jacobi polynomials. The details of the procedure considered 
here are studied extensively in references 9, 12, 13, and 14. 

Substituting equation (156) into equation (143) yields equation (157) (ref. 9): 


1_ 

n 



i s ) 


ds 




(M)<1 (157) 


The SIE (157) becomes equation (158). The terms in equation (158) are defined in equations (159), 
(160), and (161): 





ism 


in(7iX) 


+ /?„(*) 


+cX( s )\ = Qi s ) 


K fa) = J fa> 


l 


K( s ) = j {-1<S< 1) 




(158) 

(159) 

(160) 
(161) 


where (^) is the conjugate of ’^(^). 

In equation (158), c is the conjugate of c . So that the functional equation (158) can be reduced to an 
infinite system of algebraic equations with the unknown coefficients c n , both sides of the equation are 

multiplied by w(-'k,-$,s)Pf n ~ X ~^\s) with (m = 0, 1, 2, ...) and are integrated over (-1, 1). The 
orthogonality relation of the Jacobi polynomials in equation (162) is also used: 




0, n,m = 0,1,2,... 

e p„p, = 2 X+ P +1 r(m + X + l)r(m + p + l) 

2m + a + (3 + l m!r(m + ?i + P + l) 


n = in 


(162) 


Since the index of the singularity is K = 1 for m = 0, the right side of equation (162) becomes equation 
(163), with the Jacobi polynomial defined in equation (164): 


= J w{s)ds = 


2^ +p+1 r(x+i)r((3+i) 
r(x + p + 2 ) 


P 0 M) (s) = 1.0 


(163) 

(164) 
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Truncating the series of equation (158) for the first N terms yields equation (165): 

— • \ G '» + X [d mi c„ + KJn ) = T ,n {m = 0,l,..-,N) (165) 

sm^7TAj n=0 

Terms in equation (165) are defined in equations (166) to (172): 



^m + 1 CR m + 1 T i Cl m + ] 

(166) 


c n = C R n + iCI n 

(167) 


Cyi ~ CR„ ~ lCI n 

(168) 

^nm 


(169) 

1 

7* _ I 

I 


(170) 




T m = \ 

-l 

vv(-^,-(3, .v) = (1- 1 s , ) _ a '(1 + 5') _ ^ = w -1 ^) 
The matrix form of equation (165) is written in equation (173): 



Ctf] 

\ a ] 


L h(N+l),2(N+2) 

[ci\ 


2(W+2) 



where 

CR real part of c 

Cl imaginary part of c 

TR real part of the right side 
TI imaginary part of the right side 


(171) 

(172) 

(173) 


The number of equations in equation (173) is 2 (N + 1). There are 2 (N + 2) unknown constants: CRq, 

CR ] , CR n+1 and CI 0 , CI V and CI N+l . Therefore, the continuity condition is the additional equation used 
to provide the unique solution given in equation (174): 

i 

J <\>(s)ds = 0 (174) 
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Substituting equation (156) into equation (174) and using the orthogonality conditions in equation (162) 
yields equation (175): 


(ci? + iC 7 )e 0 (A,,p) = 0 

Adding an additional equation (eq. (175)) into equation (173) yields equation (176): 


( 175 ) 


[A] 


2(N+2),2(N+2) 


N 


[7-q 

(c/J 

1 | 
1 2(N+2) 

177) 


(176) 


2(N+2) 


Any standard method of solution of the system of algebraic equations can be used to determine the 
unknown c in terms of CR and Cl, as in equation (177): 


N 

[C/J: 


= [A] 


-1 

2(N+2),2(N+2) 


2(N+2) 


a 


(177) 


2(N+2) 


From equation (177), the unknown constants CRq, CR x , ..., CR N+l and CI 0 , CI X , ..., CI N+X , are known. 
These constants determine c 0 , c x , ..., c N+x , which are defined in equation (178): 

c 0 = Cr 0 +iCI 0 
c x = Cr x + iCI x 


C N + 1 — d^N+1 + iCI N+ \ 


(178) 


Gauss-Jacobi Quadrature Technique 

The objective of the numerical work in this technique involves the evaluation of the unknown 
constants. However, since the related integrals are of the Gauss-Jacobi type, these constants may be 
evaluated accurately without computational difficulty (ref. 9). The quadrate formula, which can be used 
for this purpose, is defined in equation (179) (ref. 15). 

1 N 

r ( 5 ) = J (1 “ S Y X (1 + 5 )” P f( S ) ds = X ^mf^m ) (179) 

-1 & = 1 


where s k and m- 1, . . ., N are the roots of the orthogonal polynomial P ^ ^ = 0, m = 0, . . ., N. The 

weights of equation (179) are defined in equation (180): 


Q =- 


2A-X-P+2 r(jv-x+i)r(jy-p+i) 


-^-(3 


r(w-x-p+i) 


(180) 
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( 181 ) 


K{s m )= + 

_1 n = 1 

where ^ n , w = 1, . . TV, are the roots of pj ^' ’^(^ n ) = 0. The weights of equation (181) are defined in 
equation (182): 


( 27V+X+P+2 r(7v+x+i)r(7v+p+i) 2 ^ +p 

(w+iMn+x+p+i) r(w+x+p+i) p^»(5„)/'£f l (5„) 


(182) 


Stress Intensity Factor 

In the immediate neighborhood of the crack tips, the form of the solution is identical to that of the 
plane strain, and the bonding stresses may be expressed as in equation (183) (ref. 5): 


p(x) + iq(x ) = Po{% (x) + iSR 2 (x)) 


— i (l + x 
x = co log 


sinx-jcosT 


(*-*■) 


(0 = — log 
2k 1 


( 1 + y*^ 

vW j 


l — x y 

(-1 < X < 1) 


(183) 


where the bounded functions c R l and 9t 2 are the intensities of stresses and are proportional to the external 
loads. 

Also from the formulation of the problem, the equation (|)(x) = p(x ) + iq(x ) can be written as in 
equation (184): 


^(x)=^c i w(x)p^' p) (x) 

1=1 

In equation (184), the fundamental function w(x) is defined in equation (185): 

w(x) = (1 - x) _(1/2)_, “(l + x) _(1/2)+! “ 

The term in equation (185) is defined in equation (186): 


CO = — log 
271 1 


1 1 + y 

z *" 

1-Y J 


Therefore, equation (186) becomes equation (187): 


vv(x) = ^1-x 2 j ^ 


1 + X 

1-x 


(184) 


(185) 


(186) 


(187) 
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Equation (187) can be written in different forms, as defined in equations (188) to (190): 


\ ico logf 

w(x) = ^^=e u- x 


VT” 


( 188 ) 


wi 


(x) = 


VTT 


(189) 


or 


vv(x) = 


cos T + ismi 


VT 


(190) 


This last form of the fundamental function (eq. (190)) is multiplied by Hi and uses Hi = The 
fundamental function in equation (190) is rewritten in equation (191): 


, , i cosT + zsmx 

H x ) = - / , 

1 A/l-X 


w(x) 


_ 1 (-sinx + zcosT) 


VT 


■x 


w(x) = i 


. sin T — l COS T 


VT 


(191) 


Substituting equations (184) and (191) into equation (183) yields equation (192): 


XKfT'W] 


n= 0 


SUIT — l COS T 


JR) 


= p 0 [9ti(x) + /9t 2 (*) 


SinT-i COST 


JR) 


(192) 


After necessary simplifications are made, equation (192) becomes equation (193): 

)[^1 M ■ + '^2 (*)] = X [ i£ n P N 3) ( X )] 


P 0 


(193) 


n= 0 


It can be shown that (ref. 16), in terms of (|)(x), the bonding stresses in the neighborhood of the crack 
tip may be expressed as in equation (194): 


°r~ il: rx =p(x) + iq{x) 


Gr-tirx=P + 


smx-zcosT 


VT 


r=r \x\ < 1 


(194) 
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Thus, the stress intensity (ref. 5) is defined as 


G r - 



sinx-zcosT 

Vl-x 2 


(195) 


where — 1 < jc < 1 . Since the left side of equations (194) and (195) are the same, the right side of equations 
(194) and (195) are equal to each other. Therefore, the SIFs are obtained as in equation (196): 

K(a)= 

k{ b ) = Po%{ b ) 

k 2 {a)= ~P 0 %{ a ) 

k 2 { b )= ~Po%{ b ) (196) 


Strain Energy Release Rate 


The SERRs for modes I and II for a crack along the interface may be expressed as in reference 5: 


r(n) _K (H 1 + K 1 H 3 )(H 3 + K 3 H 1 ) p 

^ _ 2 (x 1 (x 3 [(l + k 3 )|x x + (l + Kj )(x 3 ] 

gm=- 

n (yx.ii . 

2 Hi|a 3 [(l + K 3 )ni + (l + Ki)|a 3 ] 


G 2 ( b ) 


n (^ 1 + k 1 |x 3 )(^i 3 + k 3 }x 1 ) 

2 ja 1 )x 3 [(l + k 3 )(x, + (l + k,)|x 3 



(197) 


The total SERR is then obtained as 


dU(a ) 
3c 

3 U(b) 
3c 


G(a) = [k~(a) + iq(a)] 
G(b) = [k 2 (b) + k 2 2 (b)] 


n (|x 1 + k 1 |i 3 )(}x 3 + k 3 ^ 1 ) 

2 Ri(^ 3 [(l + )(Xi + (l + k 3 )(x 3 ] 
n ([Xt + K^X^ + K^) 

2 ja 1 ja 3 [(l + k 3 )jll 1 + (l + k 3 )h 3 ] 


(198) 


Application of the Method 

As an illustration of the use of the numerical approach developed in this report, the sensitivity of the 
SERR and SIF to coating stiffness and thickness are examined using material properties representative of 
polymer matrix composites. The method is applied to the case of a single fiber embedded in an infinite 
matrix with a radial load applied to the matrix. In the SIE, Q l and Q 2 are functions of closed displacement 
differences that are related to the applied load at infinity: 
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( 199 ) 


a=2(i-v 3 )K*-g;w 

Fo 

Q 2 =2(\-v 3 )^^g' 2 (x) 

Fo 


The displacement differences that are assumed to be known are associated with an elasticity solution. 
So that the solutions can be found, two different cylinders must be formed. The first is a hollow outer 
cylinder that represents the matrix. This cylinder has an infinite outer radius. Since the fiber and coating 
are considered to be perfectly bonded, the second cylinder is a solid cylinder formed by combining the 
inner solid cylinder (fiber) with the middle hollow cylinder (coating). There is no bond between the two 
cylinders, and external loads are applied separately to the two cylinders. 

In the example that is considered here, the external radial stress o rl (°o, x) = G 0 = 1 lb/in. 2 is applied to 
the outer (matrix) cylinder. Index r is for the radial direction, and index 1 is for the matrix. There is no 
load applied to the second cylinder. Then, from the elasticity solution of the problem, the surface 
displacements for the first cylinder (matrix) and the second cylinder (coating and fiber) are 


_2o 0 r 


2ct 0 Vi 

u, = — -x 


E i 


r3 = U x3 = 0 


( 200 ) 


The displacement differences that need to be closed at the matrix-coating (region 1 to region 3) 
interface can be written as 


-(u rl -u r3 ) = g l (x) 
-Ki -u x3 ) = g 2 (x) 


where r-r c and x e L. Therefore, rewriting equations in equation (201) yields 


&(*)=- 


2g 0 r C 

E , 


8 2 ( x ) 


20oV Ix 

E , 


( 201 ) 


( 202 ) 


As an example, symmetric loading and a single bond are considered where L = (-1, 1) and L ' is 
infinitely long. For the SIE, a zero stress state at infinity is solved; Q x and Q 2 are the only loadings that 
are considered. In the SIE, gi(v) and g 2 (x) are expressed as 


SiM = -2Pi 


&(*) = 2p 2 


Ho 


(i-v 3 )mm 

Ho 


( 1 -v 3 )HiH 3 


(203) 
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Equations (202) and (203) are used to find Pi and P 2 as in equation (204): 


_ 4(1 — v 3 )mn 3 G 0 r c 
_ 4(1-v 3 )^ 3 g 0 v 1 

Substituting equations (203) and (204) into equation (199) yields equation (205): 


( 204 ) 


aw= 0 

Qi ( x ) = P2 

8i( , )= faWw , 


(205) 


The effects of material properties and crack size on the SERR are now considered. In the SIE, the 
integrals from -1 to 1 that contain Jacobi polynomials are taken by using the Gauss-Jacobi quadrature 
technique. Providing an adequate number of points is important for convergence. As shown in figure 5, 
convergence begins at about 16 points. The material properties used here are the same as those used in 
reference 5. The material properties of the matrix (indicated by a subscript 1) are E { = 4.5xl0 5 lb/in. 2 and 
= 0.35, and the properties of the fiber (indicated by a subscript 2) are E 2 = 10 7 lb/in. 2 and v 2 = 0.2. 



Figure 5. — Number of Gauss-Jacobi quadrature points. 


NASA/TM — 2004-212328 


35 


The radius of the coating is 1, and the radius of the fiber is 0.5. In this example, the value of v 3 = 0.2 is 
not changed, but E 3 ranges from E x to E 2 and from E 2 to E = 10 8 lb/in. 2 The parametric studies are 
presented for the normalized mode-I and mode-II SIF: that is, equation (206) with c = L/2, and a 
normalized SERR; in other words, equation (209), which is derived from equations (207) and (208): 



( 206 ) 


G = (k 2 x +k 2 2 ) 


n 

2 


^[(l + ATj^ + fl + ATa)^] 


(207) 


7t Qii + KjVijXVij + g#t) 

2 X i iX i 3 [(1 + K 1 )h 3 + (1 + K 3 ) X i 1 


(208) 

(209) 


In figure 6, the SERR and SIF are plotted versus coating stiffness E 3 with a lower bound equal to the 
matrix stiffness E x and an upper bound equal to the fiber stiffness E 2 \ that is, E x - 4.5xl0 5 <E 3 <E 2 = 10 7 . 

The SERR and SIF are sensitive to variations in coating stiffness when the coating stiffness is close to the 
matrix stiffness. As the coating stiffness increases, the SERR increases and becomes less sensitive to 
variations in coating thickness. Figure 7 is similar to figure 6 except that the maximum coating stiffness 
(10 8 ) is 10 times the fiber stiffness. The SERR and SIF are nearly independent of coating stiffness when 
the coating stiffness exceeds the fiber stiffness. For polymer matrix composites, these results suggest that 
the driving force for crack growth will be sensitive to the stiffness of the interphase and that the driving 
force will increase significantly as the stiffness of the interphase increases beyond the stiffness of the 
matrix. For the hypothetical case in which interphase stiffness changes but the critical SERR for crack 
propagation remains constant, an increase in the stiffness of the interphase would result in crack 
propagation at a lower applied load. 

The effect of coating thickness on the SERR and SIF is shown in figure 8. The material properties 
used to obtain the data in figure 8 are E x - 4.5xl0 5 lb/in. 2 , v x = 0.35 for the matrix; E 2 = 10 7 lb/in. 2 , v 2 = 
0.2 for the fiber; and E 3 = 10 8 lb/in. 2 , v 3 = 0.2 for the coating. The ratio of the coating thickness to half of 
the bonded length (h 3 /c) is varied from 0.05 to 2. The normalized SERR and SIF versus the ratio h 3 lc is 
shown in figure 8. Since the SERR increases with increasing coating thickness, coating thickness could 
result in crack propagation at a lower applied load. This conclusion is valid only if the critical SERR is 
not affected by the coating thickness. In addition, this conclusion applies only for the material properties 
assumed in the analysis. In this analysis, the coating stiffness was greater than the fiber stiffness. A 
different analysis would be needed to evaluate the effect of interphase thickness in a material system in 
which the interphase stiffness is closer to the matrix stiffness. 
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Figure 6 . — Stress intensity factors (SIFs) and strain energy release rate 
(SERR) versus E 3 XIO 5 ranging from E-| = 4.5 xIO 5 to E 2 = 10 7 where 
c is half the bounded length (c = LI2). 
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Figure 7. — Stress intensity factors (SIFs) and strain energy release rate 
(SERR) versus E 3 XIO 5 ranging from E-j = 8x1 0 6 to E 3 = 10 8 . 
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Figure 8. — Stress intensity factors (SIFs) and strain energy 
release rate (SERR) versus /73/c for /?2 = 0 . 2 . 


Conclusions 

An analytical model was developed to calculate the SIF and SERR for an infinitely long straight fiber 
and an infinite matrix that were partly joined using a third material called the interphase (or coating). The 
addition of a third (interphase) region is an extension of a two-component model previously described in 
reference 5. The ability to evaluate the effect of interphase geometry and material properties makes the 
model more realistic for polymer matrix composites, which are known to have an interphase region. A 
preliminary set of parametric studies demonstrated the effect of interphase properties on the SERR. An 
advantage of this analytical method is that calculation of the SERR and SIF is straightforward for 
problems with cylindrical symmetry. Although cylindrical symmetry does not necessarily apply in a real 
composite structure, the method is still useful as a simple tool to provide insight about the effect of 
interphase properties on the crack driving force. The model can also be applied to other material systems 
such as ceramic or metal matrix composites by using the appropriate material properties. Further 
development of the method to account for different loading conditions and different crack types would 
make the method even more useful, particularly for experiments such as fiber pullout or growth of a 
debond area in a single fiber test where the SIF and SERR solutions are still lacking. 
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Appendix — Symbols 


Cl, Z?1 
^2> ^2> C 2 , D 2 
^ 3 ; ^ 3 > C 3 , Z ) 3 

C 

Cl, e 2 , e 3 

F 

Ci, C 2 , F 3 

Gi(?), G,(Z>), G 2 (a), G 2 (b) 
la- 1\, K 0 , K 1 
k* k* 

, /Vjj 

4 (a), 4 ( b ), k 2 (a), k 2 ( b ) 

L 

L' 

P^\* k ) 

Pi Q 
Qh Qi 

Re 


unknown constants for layer 1 (matrix) 
unknown constants for layer 2 (fiber) 
unknown constants for layer 3 (coating) 
half of the bounded length 

modulus of elasticity for the matrix, fiber, and coating 
Love strain function or Green’s function 

Green’s functions for the matrix, fiber, and coating, respectively 

known displacement gaps along the 1 and 2 directions 

strain energy release rates (SERRs) 

modified Bessel functions of order 0 and 1, respectively 

normalized stress intensity factors for modes I and II 

stress intensity factors (SIFs) 

length of the bonded interface 

length of the debonded interface 

Jacobi polynomial 

unknown stresses outside the crack, on the bonded interface 
functions of closed displacement differences along the radial and axial 
directions 

real part of a complex number 
stress intensities 


r 

r c 

r f 

r 

1 m 

S 

t 

1A/ 

w(x) 

X 


radial distance 
coating radius 
fiber radius 
matrix radius 
dummy variable 

dummy variable of the Fourier domain 
displacements along the r and v directions, respectively 
fundamental function 
axial distance 


a 


r c t , where r x - r 2 - r c 


T 

Y 

y* 

8 

e 

K 


gamma function 
r/, where r 2 = r 3 = r f 
material parameter 
Dirac delta function 
cylindrical coordinate 
index of the singularity 


(|X 1 ? Vi ), (|A 21 ^ 2 )? (ft 3 A 3 ) elastic constants of the matrix, fiber, and coating, respectively 
v Poisson's ratio 


\ dummy variable 

p variable (a function of rt) 

a n a Q ,a^ stresses along the r-, 0 -, and x-axis, respectively 

T rx shear stresses along the r-x direction 
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0 

V 2 

Subscripts: 

1 or m 

2 or / 

3 or c 


complex form of the unknown stresses outside the crack, on the bounded 
interface 

Laplacian operator in cylindrical coordinates 


related to the matrix 
related to the fiber 
related to the coating 


NASA/TM — 2004-212328 


40 



References 


1. Sih, G.C.; and Skudra, A.M.: Failure Mechanics of Composites. Elsevier Science 
Publishing Co., Amsterdam, 1985. 

2. Liebowitz, Harold: Fracture, An Advanced Treatise. Academic Press, New York, NY, 
1968. 

3. Yang, F.; and Pitchumani, R.: Modeling of Interphase Formation on Unsized Fibers in 
Thermosetting Composites. Proceedings of the ASME Heat Transfer Division, vol. 3, 
ASME-HTD-V ol. 366-3, 2000, pp. 329-337. 

4. Bogetti, T.A., et ah: Characterization of Nanoscale Property Variations in Polymer 
Composite Systems: 2. Numerical Modeling. CPSOA, vol. 30A, no. 1, 1999, pp. 85-94. 

5. Ozbek, T.; and Erdogan, F.: Some Elasticity Problems in Fiber- Reinforced Composites 
With Imperfect Bonds. Int. J. Engrg. Sci., vol. 7, 1969, pp. 931-946. 

6. Erdogan, F.; and Ozbek, T.: Stresses in Fiber-Reinforced Composites With Imperfect 
Bonding. J. Appl. Phys., 1969, pp. 865-869. 

7. Timoshenko, S.; and Goodier, J.N.: Theory of Elasticity. McGraw-Hill, New York, 1951. 

8. Muskhelishvili, Nikolaei Ivanovich: Singular Integral Equations; Boundary Problems of 
Function Theory and Their Application to Mathematical Physics. Noordhoff, Groningen, 
1953. 

9. Erdogan, F.; Gupta, G.D.; and Cook, T.S.: Numerical Solution of Singular Integral 
Equations. Methods of Analysis and Solutions of Crack Problems — Recent Developments 
in Fracture Mechanics. Theory and Methods of Solving Crack Problems. George C. Sih, 
ed., Noordhoff International Pub., Leyden, Netherlands, 1973, p. 368. 

10. Krenk, Steen: On Quadrature Formulas for Singular Integral Equations of the First and the 
Second Kind. Q. Appl. Math., vol. XXXIII, no. 3, 1975, pp. 225-232. 

11. Binienda, Wieslaw Kazimierz: Mixed-Mode Fracture in Fiber-Reinforced Composite 
Materials. Ph.D. Dissertation, Drexel University, 1987. 

12. Erdogan, F.: Mixed Boundary Value Problems in Mechanics. NASA CR-145042, 1975. 

13. Erdogan, F.; and Gupta, G.D.: Layered Composites With an Interface Flaw. Int. J. Solids 
Structures, vol. 7, pp. 1971, 1089-1107. 

14. Kuguoglu, Latife Hamide: Analytical Modeling of Imperfect Bond Between Coated Fibers 
and Matrix Material. Ph.D. Dissertation, University of Akron, 1999. 

15. Stroud, A.H.; and Secrest, D.: Gaussian Quadrature Formulas. Prentice-Hall, Englewood 
Cliffs, NJ, 1966. 

16. Erdogan, Fazil: Stress Distribution in a Nonhomogeneous Elastic Plane With Cracks. J. 
Appl. Mech. Trans. ASME, series E, 1963, pp. 232-236. 


NAS A/TM— 2004-212328 


41 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1. AGENCY USE ONLY {Leave blank) 

2. REPORT DATE 

February 2004 

3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

4. TITLE AND SUBTITLE 

5. FUNDING NUMBERS 


Fracture Mechanics Analysis of an Annular Crack in a Three-Concentric- 
Cylinder Composite Model 


6. AUTHOR(S) 

Latife H. Kuguoglu, Wieslaw K. Binienda, and Gary D. Roberts 


WBS-22-708-24-05 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E- 13923 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM— 2004-212328 


11. SUPPLEMENTARY NOTES 


Latife H. Kuguoglu, QSS Group, Inc., Cleveland, Ohio 44135; Wieslaw K. Binienda, University of Akron, Akron, 
Ohio 44325; and Gary D. Roberts, NASA Glenn Research Center. Responsible person, Gary D. Roberts, organization 
code 5150, 216-433-3244. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category: 64 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 301-621-0390. 


Distribution: Nonstandard 


12b. DISTRIBUTION CODE 


1 3. ABSTRACT (Maximum 200 words) 

A boundary- value problem governing a three-phase concentric-cylinder model was analytically modeled to analyze 
annular interfacial crack problems with Love’s strain functions in order to find the stress intensity factors (SIFs) and 
strain energy release rates (SERRs) at the tips of an interface crack in a nonhomogeneous medium. The complex form of 
a singular integral equation (SIE) of the second kind was formulated using Bessel’s functions in the Fourier domain, and 
the SIF and total SERR were calculated using Jacoby polynomials. For the validity of the SIF equations to be established, 
the SIE of the three-concentric-cylinder model was reduced to the SIE for a two-concentric-cylinder model, and the 
results were compared with the previous results of Erdogan. A preliminary set of parametric studies was carried out to 
show the effect of interphase properties on the SERR. The method presented here provides insight about the effect of 
interphase properties on the crack driving force. 


14. SUBJECT TERMS 

Singular integral equations; Strain energy release rate; Fiber-matrix interfaces; Fracture mechanics; 
Crack propagation; Stress intensity factors; Resin matrix composites; Fiber composites 


15. NUMBER OF PAGES 

48 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 


